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Abstract 

We examine the colour structure and charged particle yield for both the ti signal 
and the irreducible background processes contributing to e"*'e~ — > bbW^W~ 
production close to the ti threshold. The charged particle multiplicity for the 
various components of the cross section is computed as a function of several 
kinematic variables. Our study may have important implications for recently 
proposed studies of interconnection phenomena in tt production at high-energy 
e+e~ colliders. 



1. Introduction 



One of the most important physics topics at future e^e" hnear coUiders is the detailed 
study of the top quark. In particular, a unique precision in the top mass determination 
(with an accuracy of a few hundred MeV) is anticipated, see for example Refs. 0, ^ |[ . 
One of the obvious requirements for the success of such precise studies is a detailed 
understanding of the background processes, in particular of the colour structure of the 
e^e~ bbW^W~ background events. 

Recall that the dominance of the t bW~^ weak decay mode leads to a large top 
width Tt, which is approximately 1.6 GeV for ^ 180 GeV. This width is larger 
than the typical hadronic scale /i ~ 1 fm~^, and so the top decays before it has time 
to hadronize [|5|, |[. It is precisely this large decay width that makes top physics so 
unique. First, the top decay width Tf provides an infrared cut-off for the strong forces 
between the t and the i |^, ^, 0: the width acts as a physical 'smearing' [10|, and top 
production becomes a quantitative prediction of perturbative QCD, largely independent 
of non-perturbative phenomenological algorithms. Second, Tt controls the size of QCD 
interferences between radiation occurring at different stages of the ti production process 



11, O, Ol. These interferences affect the structure of the colour flows in the tt events, 



and may provide a potentially serious source of uncertainty in the reconstruction of the 
final state. 

As is discussed in detail in Ref. [|14| (see also flSl), particle production in events of 
the type 

e+e" ^tt^ bbW+W- (1) 

could depend in a non-trivial way on the kinematics of the process. The effects of ener- 
getic perturbative gluon radiation are suppressed because of the space-time separation 
between the decays of the t and i quarks. Soft gluon emission and non-perturbative 
fragmentation introduce a correlation (interconnection) between the b- and 6-jets com- 
ing from the top decays. As a result, the b and b have to 'cross-talk' in order to produce 
a final state made up of colourless hadronsQ. 

Such interconnection phenomena could affect the final state in tt events in many 



respects. In order to quantify their size, it was proposed in Ref. WM to examine the 

^Analogously to Ref. Q, throughout this paper we assume for simplicity that the W's decay 
leptonically. 
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average multiplicity of double leptonic tt decays as a function of the relative angle be- 
tween the b- and the 6-jets. Such studies could provide important information about 
the size of the interconnection-related systematic uncertainties in the top mass recon- 
struction. They are also interesting in their own right, since they potentially open 
up a new laboratory for probing non-perturbative QCD dynamics, see for example 
Refs. ig,|ll,|16|,|l3. 

The issue of background processes was not addressed in Ref. \T^\. Since the inter- 



connection effects in the distributions of the final-state particles in the process 

e+e- bbW+W- (2) 

are comparatively weak, and since the background contributions are quite different 
in various kinematical regions, the detailed properties of the colour flow structure of 
the background events need special detailed consideration. Here we present a com- 
plete analysis of the full process (0) by studying (i) the detailed colour structure of ti 
production and decay, and (ii) the colour structure of the irreducible backgrounds in 
bbW~^W~ production. 

This paper is organised as follows. In Section 2 we give a detailed description of the 
structure of the matrix element for the reaction e+e^ bbW^W~. Numerical results 
are given in Section 3 and the conclusions are presented in Section 4. 

2. Colour Structure of the Matrix Element 



The matrix element used in this study is the one already introduced in Refs. [0, |19 
where further details (method of calculation, parameter values, etc.) can be found. 
Only the value of the Higgs mass Mfj has been changed in the present work. Since 
in the present context we are not interested in effects due to virtual Higgs states, we 
deliberately increased Mh to a value greater than the collision energy, viz. Mh = 
600 GeV for y/s < 500 GeV. We use the value nit = 175 GeV for the top quark mass. 



as suggested by recent data |20 



With the assumption that the Higgs mass is always larger than the collider CM 
energy, the dominant contributions to the modulus squared of the scattering amplitude 



are (see Ref. [19 



(i) double t-resonance (i.e. tt bbW^W , Fig. la). 
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(ii) single t-resonance (i.e. t — > bW^ + charged conjugate, Fig. lb), 

(iii) Z°-resonance (i.e. —* bb, Fig. Ic), 

(iv) non-resonant (NR) (Fig. Id together with those in Fig. Ic involving the splitting 
process 7* — >• bb). 

If we denote by Mj the sum of the diagrams in the jth channel, and neglect Higgs 
contributions, we have 

4 

Mtot = J2Mj, (3) 
i=i 

where Mtot is the total Feynman amplitude. In squaring Eq. (^, we define the combi- 
nations0 

Ml = |Mi|2, (4) 
M2 = |M2|2 + 2Re[MiM*] (5) 

^3 = |M3|^ (6) 

A^4 = |M4|^ + all remaining interference terms, (7) 



such that 
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We associate these four combinations with the contributions (i) - (iv) defined above. 

When discussing the structure of the particle flows corresponding to the background 
processes (ii), (iii) and (iv), it is convenient to apply the standard 'parton shower plus 
fragmentation' picture for e+e^ ^ 7*, Z° ^ qq, various aspects of which are now well 
understood. Analogously to Ref. the colour flow distribution can easily be de- 
scribed using the language of QCD antennae/dipoles (for details see [|^). Thus, the 
radiation pattern corresponding to the diagrams of Fig. Ic and Id is given by the bb 
antenna which describes the production of hadrons in the process 7*, 2'^ — > bb. Because 
of the large amount of data from PETRA/PEP and LEPl/SLC, the properties of 6-jets 
in this process are now well understood. A specific role is played here by kinematical 

^Note that only the sum of these is positive definite. In particular, the interference terms in Aii 
can cause the corresponding contribution to the cross section to be zero or negative, see Figs. 3 and 4 
below. 
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configurations in which the b- and 6-jets come from the decay of an on-mass-shell 
boson. For these configurations the average particle multiplicity is given by the multi- 
plicity in bb, as measured for example at LEPl. For the non-resonant processes 
in Fig. 1 the average multiplicity of the final-state particles depends strongly on the rel- 
ative orientation of the b- and 6-jets. In contrast, for the double t-resonance (Fig. la) 
and single t-resonance (Fig. lb) mechanisms the angular dependence is comparatively 
weak (but nonetheless readily visible) because of the suppressed role of the bb antennae 



0. 



It is worthwhile to mention that not far from the ti threshold the final-state struc- 
ture corresponding to double top-resonant production is dominated by two essentially 
azimuthally symmetric b- and 6-jets, while the background 7*, 2'° bb processes are 
described by the bb antenna patterns with all the appropriate azimuthal asymmetry 



effects (see Ref. [^). Moreover, the 7* — > bb contribution enhances the role of low 
multiplicity events. Note also that with increasing collision energy the b- and 6-jets 
coming from tt production tend to be more back-to-back, while those resulting from 
7*, bb splitting tend to become more coUinear. 

3. Results 

Our numerical results are presented in Figs. 2-8. In all of them we distinguish between 
the e~^e~ —>■ bbW^W~ components (i) — (iv) discussed in the previous section and in 
Refs. lT|,|g. 



In Fig. 2 we show the dependence of the various contributions on the collider centre- 
of-mass (CM) energy, in the range 200 GeV ~ ^/s ~ 500 GeV. The behaviour of the 
various curves can be readily understood in terms of the diagrams contributing to 
the J^] amplitudes, see Fig. 1. The ti bbW~^W~ curve corresponds to tt-pair 
production followed by the decays t bW^ and t — > bW~ , both above (i/i > 350 
GeV) and below threshold {^/s < 350 GeV). In contrast to the naive case of the 'narrow 
width approximation' (see Ref. [0), in which this process is computed in terms of a 
productionx decay reaction a{e^e~ — > ti) x [BR{t — > 6W^)]^ with zero top width in 
the production process (top quarks produced on-shell) and independent semileptonic 
top decays, here both t-quark finite width effects and energy/spin correlations between 
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the two top decays are included. The consequence^ of this is clearly visible in the 
dependence of the cross section on y/s. In particular, there is a significant 'shoulder' 
which extends far below the nominal threshold. 

This behaviour naturally also affects the other e~^e~ — > hhW^W~ contributions that 
include interferences with the double top resonance channel, such as t — > hW^ and NR. 
This is clearly visible in both of these as the little 'bump' around ^/s = 2mt, which 
slightly decreases the cross section in the first case and enhances it in the second. It 
is important to stress that these curves, and therefore their interplay, have no separate 
physical meaning, since the A1?'s are not separately gauge invariant, but they provide 
a useful way of looking inside the process and distinguishing between the different 
interactions. The final step must always be to sum together all the various components 
M-f and consider the total rates. 

On their own, the curves in Fig. 2 corresponding to t ^ hW^ and NR describe, 
respectively, the production of a single top quark through a highly virtual W^* decaying 
to th pairs with t hW , and dominantly (as will be apparent from Figs. 3-4, see 
below) of a 66-pair from 7* bb. In the first case the main contribution to single top 
production comes from diagrams in which the decaying is produced via t-channel 
exchange (i.e. the second diagram of Fig. lb) rather than via s-channel exchange 
(i.e. the first and third diagrams of Fig. lb), since the corresponding curve in Fig. 2 
does not show any decrease with ^/s above the nominal threshold at rrit + My/ + 
rrib ^ 260 GeV. In the second case, except for the range 270 GeV ~ y/s ~ 360 GeV 
where the interferences between non-resonant and resonant diagrams are important, the 
shape is almost flat, reflecting the dominance of the t-channel contribution (the first, 
second and fifth diagrams of Fig. Ic, with a 7 propagator) to 'j*W~^W^ production 
followed by the 7* bb splitting. Finally, the curve corresponding to the 
bb background illustrates the characteristic features of the production of three vector 



Other effects are discussed in Rcf. [18 . Note that analogously to Ref. we do not address 
here the issue of the QCD Coulomb attraction effects in process (|l|). These effects are especially 
important near the ti threshold where they substantially enhance the production amplitude. The 
Coulomb physics could be incorporated using the technique of non-rclativistic Green's functions for 
the tt system |Q, ^ which has now become quite routine, see Refs. |l|, ||, ||, Q and references therein. 
The background amplitudes which are our main concern in this paper are only weakly affected by 
QCD final-state interactions. We should also mention that in the present study we do not take into 
account QED Initial State Radiation (ISR), but these effects would be straightforward to include. 
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bosons Z^W~^W~ , with the Z° decaying to the fe6-pair, with a threshold at 2Mw+Mz ~ 
250 GeV and a dominant t-channel exchange structure (the first, second and fifth 
diagrams of Fig. Ic, with a Z°-propagator). 

As already mentioned, we concentrate on the case of leptonic decays of both the 
W^^s. This implies that we have final states with two neutrinos and therefore with 
missing energy and transverse momentum: bbW~^W~ — > 66(£z/^)(£V£/) or more generally 
jjW^W" jj{^vi){^'ve.')i depending on whether or not one can exploit the possibility 
of using 6-tagging This missing energy and transverse momentum prevents us 

from studying the kinematics of e^e'^ bbW^W~ events by looking at variables which 
require the reconstruction of the VT^'s through their decay products, for example the 
invariant masses Mi^w or energies -EfeVK- These quantities would of course be extremely 
useful in distinguishing the ti signal from the other e^e~ bbW^W~ background 
contributions [0, 0. Thus in the figures which follow we restrict our attention to 
variables involving the 66-system only. 

In Fig. 3 we show the dependence of the four components (i)-(iv) on the cosine of the 
angle between the b and 6, at ^/s = 360 GeV. The convenience of this variable has been 



explained above (see also Ref. IQ)- For double and single top production the angular 
dependence is weak, refiecting the absence of any strong correlation between the top and 
anti-top decay products near threshold. In contrast, the bb and NR components 

exhibit a more pronounced structure. For the former, we note that there is a significant 
boost to the at this centre-of-mass energy (note that ^/s = 360 GeV is 110 GeV 
above the Z^W^W~ threshold), and so the peak at cos(66) = —1 corresponding to 
back-to-back production in the Z^ rest frame is significantly smeared out. For the NR 
component, the main feature is the 7* — » bb contribution which prefers small invariant 
masses and therefore has cos(66) ~ 1. 

In Figs. 4 and 5 we show the dependence of the four components on the invariant 
mass of the 66-system and on the energy of the b{b) respectively. As expected, for the 
former we can clearly recognize both the M^^i Mz peak in the Z^ bb component 
and the 7* bb contribution at small in NR. Neither the double top nor the single 
top channel show any particularly distinctive features. Fig. 5 indicates the average 
energy of the 6's coming from the various subprocesses: approximately 70 GeV for 
tt bbW+W-, 80 GeV for tbW+ and 60 GeV for Z° bb. For the NR component 
there is a preference for small b energies. 
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We next adopt a semi-qualitative procedure for evaluating the average charged 
multiplicities corresponding to the different contributions to the e+e^ bbW^W' 
events not far from the ti threshold. This procedure is based on the antenna pattern 
analysis [|ll|, |l^, |l3l of soft gluon radiation in top quark production events within the 



framework of the MLLA-LPHD picture of multiple hadron production 

Our starting point is the multiplicity Njji{W) in e~^e~ —>■ bb production at CM energy 
W = ^/s = 2Ei,, where is the energy of the &-jet in the CM system. According to 



the MLLA-LPHD scenario (see Ref. [23| for details) A*"^^ can be expressed in terms of 
the multiplicity Ngg in the light-quark production process e'^e^ — » qq as 

N,i{W) = N,g{W) + St, (9) 

where the difference 5^ is ly-independent. This prediction of perturbative QCD is in 
a good agreement with existing experimental data |2^, from which one obtains 6b — 3. 
For Nqq{W) one can use the MLLA-inspired simple analytical formula proposed in 
Ref. |^5|. As a result, the charged particle multiplicity in bb events can be readily 
evaluated asQ 



Ar|°'^^(E,(-,), cos(66)) ^ ATf '^^(M,5) « 5.85 + 0.125 exp[2.317v/log(M,5/GeV)], (10) 

withg M?; = 2{EtEi - \pM cos{bb) + ml) and = E^^^ - ml 

Motivated by the results of Refs. ITT], |1^, we can approximate the average mul- 
tiplicities corresponding to the t-resonant production processes ((i) and (ii)) by 

Nll'\Et(i),cosibb)) ^ Nll'\Eb,Ej,) ^ ^ATf '^^(M,-, ^ 2E,) + ^Arf '^^(M,^ ^ 2^-,), 

(11) 

where Eb and Ei are the actual CM energies of the 6- and 6-jets, respectively. This 

result accounts for the contributions of the antennae of the type tb, ib or tb, ib (note 

that interconnection effects are not included here). 

■*0f course one could alternatively apply all the powerful machinery of the successful Monte Carlo 

algorithms (JETSET, HERWIG, ARIADNE). However for the purposes of this paper we can safely 

restrict ourselves to the estimates based on Eq. (p^). It is remarkable that in a wide energy domain the 

predictions of Eq. (|l0| ) are in a very good agreement with the results of the JETSET parton shower 

Monte Carlo program (We are grateful to T. Sjostrand for providing us with the numerical results 

of JETSET for N^i{W)). 

^The superscripts tt, t, Z'^ and NR in Eqs. ( |l0| , |rTl ) refer to the production mechanisms (i)-(iv) of 

the 6fe-pair in e+e~ bbW^W~ events. 
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Using the above formulae (i.e. Eq. ( p!OD for the multiphcity for the bb and 

the non-resonant contributions ((iii) and (iv)) and Eq. (|IT]) for the double and single 
t-resonant contributions ((i) and (ii)), we can study the quantities 

(ATi X _ IdE,^i^T{E,^j,),cos{bb)) 

^ JdE,^,^J'^{E,^,),cos{bb)y ^ ' 

,J.r^^ _ /c?cos(&&)P(Efc(^),cos(&&)) 

^ /rfcos(66) J^(E,(5),cos(66))' ^ ' 

at fixed CM energy (i/i = 360 GeV), and 

. ^ JdEi,(i^dcos{bb) cos(6&)) 
^ /dE,(5)rfcos(66):r^(E,(5),cos(66))' ^ ^ 

as a function of the CM energy in the range 240-390 GeV, with i = tt, t, and NR. 
In Eqs. (p^-(p^ J* and J"^ are given by 



(15) 



-6(6) 

and 



:r^(i5;,(5), cos(66)) = — ^— -=r, (16) 

^ ' dEfj(i-^d cos{bb) 

where a* is the total 'cross section' of the i-th component. 

The above quantities can be interpreted as the average number of charged particles 
per unit interval of the cosine of the angle between the b- and 6-jets (Eq. (|T^)) and 
per unit interval of the energy of the 6(6)-jets (Eq. (p!3|)), and as the total number of 
charged particles for all cos (66) angles and E{b) energies (Eq. (p!4D). These are defined 
for each of the four components (i)-(iv) separately. The distributions for {N^qs), {N^j) 
and (A^tot) plotted in Figs. 6, 7 and 8, respectively^. 

In Fig. 6 the multiplicity arising from the double and single t-resonant contributions 

is essentially independent of the angle between the 6-jets. This follows from the fact 

that (i) the multiplicity here (Eq. (pT])) does not depend explicitly on the angle, and 

(ii) the distribution of the energy of the 6-quarks (Fig. 4) is rather sharply peaked 

in both cases around 60 — 80 GeV. The multiplicity from the Z^ bb contribution 

is also approximately independent of the angle (reflecting the fact that M^g ^ Mz is 

^Note that in Figs. 6-8 a cut on the energy of the b- and &-jets has been apphed, i.e. E{b, b) > 20 
GeV. 
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constant) apart from near cos{bb) = 1 where M^j < Mz- In this region, however, the 
corresponding cross section is very small, see Fig. 3. The non-resonant contribution 
decreases steadily as cos(66) increases as the 7* bb contribution (which prefers small 
Mf^l and hence smaller multiplicity) becomes more important. 

In Fig. 7 the multiplicity arising from the double and single t-resonant contribu- 
tions rises smoothly with increasing Eh, as expected from Eq. ([Tl|). The slight fall off 
at very large Ei, (where the cross section is very small, see Fig. 5) is caused by asym- 
metric configurations where Ei, ^ E^ or vice versa. The multiplicity from the bb 
contribution is essentially constant over the whole energy range and the non-resonant 
contribution multiplicity increases quite sharply with Eb, reflecting the underlying de- 
pendence on M^i- 

The collider energy dependence of the multiplicities arising from the various con- 
tributions is shown in Fig. 8. The double and single t-resonant contributions increase 
approximately linearly with ^/s, reflecting the underlying increase in the energy of the 
6-quarks. The bb contribution is approximately independent of ^/s above the 

threshold at 2Mw+Mz ~ 250 GeV for producing on-shell Z° bosons. The non-resonant 
multiplicity is somewhat smaller, reflecting the on-average smaller values of M^i, and 
increases slowly with -y/i. Also shown (continuous line) in Fig. 8 is the sum over the 
four components, obtained using the formula 

(iVtoT) = ^^pT^, (17) 

where i refers as usual to the contributions tt — > bbW~^W~, t bW^ + charged 
conjugate, Z^ bb and non-resonant. This total multiplicity simply interpolates 
between the non-resonant contribution, which dominates the total cross section for 
~ 250 GeV, and the double-resonant tt contribution, which dominates the total 
cross section from just below the tt threshold, see Fig. 2. 



4. Conclusions 

A high-energy e^e" collider would provide the ideal laboratory for precision top-quark 
physics. Whether one considers an accurate measurement of the top-quark mass or 
a detailed study of the non-perturbative QCD dynamics and interconnection effects 
associated with tt production, one certainly needs a more detailed understanding of 
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the colour structure of the background events in order to separate their influence from 
the subtle interference phenomena in the top signal. 

In this paper we have examined the colour structure and charged particle yield in 
both the tt signal and background processes contributing to e+e" — > bbW'^W~ produc- 
tion not far from the tt threshold. Since the relative contribution of the background 
processes to the total particle yield in e'^e~ — > bbW'^W~ events is relatively small (see 
Fig. 2), we have been able to use an approximate procedure for evaluating the charged 
particle multiplicities. 

Figures 6-8 contain the main results of our study. They show the multiplicities 
associated with the various resonant and non-resonant components of the cross section, 
as functions of the energies of the b- and 6-jets, the angle between them, and the 
overall collider energy. The colour structure of the various components is very different 
(see Fig. 1) and this is reflected in the multiplicity distributions. It is important, 
therefore, that the background contributions are carefully taken into account in studies 
of interconnection phenomena in tt production. 
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Figure Captions 



[1] Typical Feynman diagrams describing at tree-level the four e+e" bbW'^W~ 
components: (i) double i-resonance (a), (ii) single i-resonance (b), (iii) Z^- 
resonance (c, with Z°-propagators) , and (iv) no resonance (c, with 7-propagators, 
and d). Internal wavy lines represent 7, or W^, as appropriate. Permutations 
of (1^+, 7*/Z°*) and (7*,^°*) along the fermion hues, exchanges between 
<-> W~ in graphs involving three-vector-vertices, as well as the charged 
conjugate diagrams of (b) are not shown. 

[2] Dependence of the individual "total" cross sections o"j on the collider energy y/s 
for the four e+e~ bbW~^W~ components: (i) double i-resonancc (continuous 
line), (ii) single i-resonance (dashed line), (iii) Z°-resonance (dotted line), and 
(iv) no resonance (chain-dashed line), for rrit — 17b GeV and Mh — 600 GeV. 

[3] Differential distributions dai/dcos{bb) in the cosine of the angle between the b 
and b quarks for the four e~^e~ — > bbW^W^ components: (i) double t-resonance 
(continuous line), (ii) single t-resonance (dashed line), (iii) Z°-rcsonance (dotted 
line), and (iv) no resonance (chain-dashed line), at ^/s = 360 GeV, for nit = 175 
GeV and Mh = 600 GeV. 

[4] Differential distributions dai/dM{bb) in the invariant mass of the 66-pair for the 
four e+e~ — > bbW'^W~ components: (i) double t-resonance (continuous hne), 
(ii) single t-resonance (dashed hne), (iii) Z°-resonance (dotted hne), and (iv) 
no resonance (chain-dashed line), at y/s — 360 GeV, for rrit — 175 GeV and 
Mh = 600 GeV. 

[5] Differential distributions dai/dE{b) in the energy of the b(b) quark for the four 
e+e~ — > bbW'^W~ components: (i) double t-resonance (continuous hne), (ii) 
single t-resonance (dashed line), (iii) Z°-resonance (dotted line), and (iv) no 
resonance (chain-dashed hne), at = 360 GeV, for rrit — 175 GeV and Mh — 
600 GeV. 

[6] Dependence of (Ncos), defined in the text, on the cosine of the angle between the 
b and b quarks for the four e^e~ bbW~^W~ components: (i) double t-resonance 
(continuous line), (ii) single t-resonance (dashed line), (iii) Z°-resonance (dotted 
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line), and (iv) no resonance (chain-dashed hne), at ^/s = 360 GeV, for nit = 175 
GeV and Mh = 600 GeV, after imposing the cut E{b, b) > 20 GeV on the energy 
of the 6's. 

[7] Dependence of (Ne), defined in the text, on the energy of the b(b) for the four 
e~^e~ —>■ bbW^W~ components: (i) double t-resonance (continuous line), (ii) 
single t-resonance (dashed line), (iii) Z°-resonance (dotted line), and (iv) no 
resonance (chain-dashed line), at ^/s = 360 GeV, for rrit = 175 GeV and Mfj = 
600 GeV, after the cut E{b,b) > 20 GeV on the energy of the 6's. 

[8] Dependence of (A^tot), defined in the text, on the collider energy i/i for the four 
e~^e~ — > bbW~^W~ components: (i) double t-resonance (continuous line), (ii) 
single t-resonance (dashed line), (iii) Z°-resonance (dotted line), and (iv) no 
resonance (chain-dashed line), and for their sum (see Eq. (p!7|)), for nit = 175 
GeV and Mh = 600 GeV, after imposing the cut E{b, b) > 20 GeV on the energy 
of the 6's. 
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